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Abstract. - We show that a mesoscopic system such as Feynman's ratchet may operate as a 
heat pump, and clarify a underlying physical picture. We consider a system of a particle moving 
along an asymmetric periodic structure. When put into a contact with two distinct heat baths of 
equal temperature, the system transfers heat between two baths as the particle is dragged. We 
examine Onsager relation for the heat flow and the particle flow, and show that the reciprocity 
coefficient is a product of the characteristic heat and the diffusion constant of the particle. The 
characteristic heat is the heat transfer between the baths associated with a barrier-overcoming 
process. Because of the correlation between the heat flow and the particle flow, the system 
can work as a heat pump when the particle is dragged. This pump is particularly effective at 
molecular scales where the energy barrier is of the order of the thermal energy. 



Introduction. - In the last decade, physics and technology have moved more and more 
toward the manipulation of mesoscopic objects. Techniques of micro manipulation of molecules 
are expanding to open new possibilities to future technology. Although to control heat flow 
in mesoscopic scales may be among crucial issues, its study is still in a preliminary stage [1]. 
Heat pumps working in nano-scale devices [2, 3] or mesoscopic machines [4] are studied but 
the number of cases are still limited. In this Letter, we study a class of mesoscopic systems 
working as a heat pump and present a physical picture explaining how they work. We start 
from a numerical study adopting one specific model, and demonstrate that a characteristic 
heat in a barrier-overcoming process determines the ability of the heat pump and Onsager's 
reciprocity coefficients. Then, we give a theoretical derivation of the reciprocity coefficients 
and show the generality of the results in a class of Feynman's ratchet [5]. We present a simple 
and clear picture for the mechanism of heat pumps, which is expected to be useful in designing 
various mesoscopic heat pumps. 

Model. - We study a class of systems including Feynman's ratchet in which a particle 
moves along a periodic structure. The particle interacts with the periodic structure in an 
asymmetric manner. The periodic structure and the particle are attached to distinct heat 
baths Bi and B2, whose temperatures are Ti and T2, respectively. We start from a specific 
but typical model shown in Fig. ^ [6], which is suitable for numerical simulation. In this 
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Fig. 1 - (a) Profile of tfie model for numerical simulation. We fix 7 = 1, Nc — 7 with a periodic 
boundary condition, and Kc = 0.5 except for Fig. 2] The depth of Vas{xp — Xi) is 0.2. (b) The flow 
of energy in the system. Ji and J2 are mean heat flows from heat baths Bi and B2 to the system, 
and W is the mechanical work per unit time provided by the external force. 



model, a one-dimensional harmonic lattice corresponds to the periodic structure. The time 
evolution of the system is described by a set of Langevin equations, 

X. = -Ji, + m - ^^^"^ , = -7ip + CpW - 1^ + /, (1) 

where Xp and Xi are the positions of the particle and the i-th lattice site (i — I, ■ ■ ■ , Nc). An 
external force / is applied to the particle. 7 is the friction constant and £,a{t) is the random 
force satisfying {^a{t)) = and {^a{t)^i3{t')) — Sa^pS{t — t'). The interaction between the 
particle and the lattice sites is described by the asymmetric potential Vas = J2i'^asixp — Xi) 
with Was as in Fig. ^ M = X]i^o(a^i — 2:^0) + Vcixi+i — Xi) is the harmonic potential for 
the lattice, where Vo{x) = Kcx'^/2 is the on-site potential and Vc{x) — Kc{x — /)^/4 is the 
interaction between neighbors. / is the lattice interval, Xio = il for i-th site, and is a 
parameter for stiffness. When Kc is sufficiently large, the lattice does not fluctuate and the 
potential for the particle reduces to a one-dimensional sawtooth-shaped potential. 

On this periodic structure, the particle performs stepwise motions with a typical step size 
At sufficiently low temperatures, each stepwise motion occurs stochastically with a dwell time 
long enough for the particle to lose its memory. In equilibrium, the probabilities of rightward 
and leftward steps are of course identical. In the following, we concentrate on nonequilibrium 
steady states maintained by an external force and/or a temperature difference. 

Heat flow induced by particle motion and Onsager relation. - Temporal heat flows from 
the heat baths to the system are defined for each trajectory by using the stochastic energetics 
[7,8] as 

•^1 W = + V^^iUi)] ° Mt) = h7ip + v/2^CpW] ° ip, (2) 

i 

where Ji{t) is the heat flow from Bi to the system and J2{t) is that from B2 to the system. 
From Eq.lPJ, we define temporal heat flow between the two heat baths as 

Ut)^{Mt)-Mt))/2. (3) 

In the steady states, mean heat flows Ji = (Ji(i)), J2 = {J2{t)) and Jq = {Jq{t)) can be 
defined, where (•) denotes average over time and ensemble. 
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Fig. 2 - (a) Mean heat flows Ji and J2 as a function of f at Ti — T2 — 0.05. For small /, either Ji 
or J2 is positive, which indicates that the system absorbs energy from one of the baths, (b) Jq vs / 
for Ti = r2 = 0.05 and 0.033, where Jq oc /. 



When a steady external force / is applied to the particle, energy corresponding to the 
mechanical work by / is provided to the system. Its mean rate is 

W = (xpf) = Jp/, (4) 

where Jp is the mean particle flow defined as Jp = (ip). In the steady states, the energy of 
the system remains constant on average. This constraint leads to the balance of the energy 
flows W, Jl and J2 as in Fig. ^h), and imphes W = —{Ji + J2)- 

Let Ti —T2. When the particle is dragged, the particle and the sites are expected to warm 
up due to the friction. Then one expects that heat transfer from the hot components to their 
surroundings takes place, and that Ji < and J2 < 0. On the contrary to this expectation, 
Jl or J2 can be positive for small values of / in Fig. EJa) (^). The system then works as a 
heat pump which absorbs energy from one of the heat baths and dissipates it to the other. 

Fig. |2Ib) shows that the heat flow Jq depends linearly on / in a wide range, including the 
region with Ji J2 > (see Fig. ^a.) ). The particle flow Jp also shows a linear dependence 
on / in a comparable range. Since the range in which the system works as a heat pump is 
included in this linear range, we are led to study Onsager relation [9] for Jp and Jq. Define T 
and AT by Ti = T - AT/2 and T2 = T + AT/2. Thermodynamic forces for the flows Jp and 
Jq are //T and AT/T^, respectively. Onsager relations are formally written as, 

_ / AT _ / AT 

= ^pp^ + -^pq"^^' "^q ~ -^qp^ + -^qq"^^- (5) 

The coefficients are written as 

Lpp^Dp, L^^^T^X, L^^ = Lp^ = q-'iDp/l, (6) 

where Dp is the diffusion constant for the particle in equilibrium, A = dJq/dAT is the heat 
conductivity of the system and q°'^ is the characteristic heat which will be discussed in the 
next section. To derive the expression for Tpp, we define the mobility as /i = dJp/df. Then 

(^)If the system is not in overdamped limit, the mechanical force can make some part of the system cool 
down. When J2 > 0, kinetic temperature of the particle is lower than T2 and, when Ji > 0, that for the lattice 
sites lower than Ti. 
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Fig. 3 - Schematics of energy flow when the system overcomes an energy barrier in equilibrium. 
The bright and dark arrows indicate the heat transfer from the heat bath Bi and B2, respectively. 
In the rightward step, Qi is transfered from Bi during the uphill and —Qi during the downhill. 
Because of the time-reversal symmetry in equilibrium states, Q'l is transfered from Bi during the 
uphill and ~Qi during the downhill in the leftward step. Similar relations hold for B2 and the 
relation Qi + Q2 ~ Q'l + Q'2 holds. Because the system has an asymmetry, Qi 7^ Q'l and Q2 7^ Q2 in 
general. 



Einstein's relation /i = Dp/T implies Lpp = /^T = Dp. The expression for Lqq follows from 
the definition. The expressions for the off-diagonal coefficients Lqp and Lpq will be derived in 
the following (see and ifTUI) ). 

Characteristic heat associated with barrier- overcoming process. - In this section, we 
investigate a characteristic heat associated with a barrier-overcoming process in the vicinity 
of the equilibrium states. It enables us to determine the non-diagonal coefficients Lqp and 
Lpq. First, let us consider equilibrium cases. When the system exhibits a stepwise motion, it 
should overcome an energy barrier. To reach the top of the barrier the system absorbs some 
amount of heat from the two heat baths and dissipates the same amount of heat to descend 
from the top (Fig. In equilibrium, of course, no net heat transfer occurs between the baths 
on average. But if we take the ensemble averages of the heat transfer for the individual steps 
classifying the rightward and the leftward steps, nonvanishing heat transfer between the two 
baths can survive. We find that a single rightward step carries heat q^, and a leftward step 
carries on average. In equilibrium, the rightward and the leftward stepwise motions occur 
with the equal probability and the relation = — gL = q°'^ holds, which consistently leads to 
vanishing of net heat transfer between the two heat baths. In general, nonvanishing values of 
g'"^ are obtained [10], because it comes from left-right asymmetry of the system. 

Next, consider the case where a small external force / is applied under the iso-thermal 
condition AT — 0. From the measurement of Jp and Jq, we can define a characteristic heat 
Qc = Jq/{Jp/l) in the hnear region of Jp and Jq. Because Jp/l is the net number of directed 
steps per unit time, qc is conjectured to be a heat transfered per single directed step. By 
determining the value of qc and g*"^ separately from numerical simulation, we have observed 
qc ~ for sufficiently small / (see Fig. Ela)). This implies that the characteristic heat qc 
associated with the stepwise motion in the linear region is the characteristic heat q""^ at the 
equilibrium states. Then, we obtain 



J. = ^f-'-j^^, (7) 



by substituting Jp — /// into Jq/(Jp/Z) = q"^^. Eq. holds for a wide range of the model 
parameter Kc as is examined in Fig. Elb), where Xq = Jq/ f- Thus the non-diagonal coefficient 
Lqp is obtained as q'^'^Dp/l. In addition, another non-diagonal coefficient Lpq briirgs the same 



Naoko Nakagawa and Teruhisa S. Komatsu 



A HEAT PUMP AT A MOLECULAR SCALE CONTROLLED BY A MECHANICAL FORCE f 




Fig. 4 - (a) Replot of the data in Fig|3;b). Both data for T = 0.05 and 0.033 are scaled to the same 
form as Jq ~ q'^'^Jp/l, which corresponds to Qc — q'^'^. Jq and Jp are calculated for / 7^ while q^'^ is 
evaluated for / — 0. (b) The agreement of Xq with fiq'"^/l over various Kc. T = 0.05. Xq = Jq/f 8'iid 
^ = Jp/ f for small values of / with AT — 0. q'^'^ is evaluated for / = and AT = 0. At the same 
time Txp is plotted, where Xp = Jp/AT for small values of AT with / = 0. The observed equality 
Xq = Txp indicates that the reciprocity relation Lqp = Lpq holds. 



form with iqp via the agreement Xq = T'Xp Fig. Elb), where Xp = Jp/^T for small AT 
with f — 0. We can rewrite the reciprocity coefficients as Lqp = ipq — q'^'^n'^'^l using the 
rate constant k''^ for the barrier-overcoming process of the stepwise motion in equifibrium, 
where Dp—K°'^P. This form clearly demonstrates that the characteristic heat brought by the 
barrier-overcoming process results in the reciprocity coefficients. 

Generality of the form for reciprocity coefficients. - In this section, we show that the 
expressions of the reciprocity coefficients are valid in a class of Feynman's ratchet attached 
to two heat baths. In this class of systems, we can generally define particle flow Jp and heat 
flow Jq between the two baths and those temporal quantities ip{t) and Jq{t). The reciprocity 
coefficients are written with the time correlation function in equilibrium states [11] as 

Lpq = Lqp = lim — / dt / dt' {Xp{t)Jq{t + t')) (8) 

where (•) denotes ensemble average. Nonvanishing contribution in the r.h.s. of Eq. JSJ comes 
from the stepwise motion of the particle. Supposing that the correlation decays within a finite 
time and that each stepwise motion is sufficiently separated in time, Eq. (jS)) is rewritten as a 
summation of the contribution from each stepwise motion at t — ti^ as 

E f^^^dt{xp(t)fdt'jq{t + t')), (9) 

where is is the index of the stepwise motion and r is a time sufficiently longer than the 
correlation time. Because the integration cji^ = dt'Jq{t + t') has no expficit dependence on 
t, the integration of Xp with t can be done separately, which gives +1 and —I for rightward 
and leftward stepwise motion, respectively. Classifying the stepwise motion indexed by ig into 
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Fig. 5 - The dark shaded area shows the range where the system functions as a heat pump when 
T — 0.05. Numerically identified boundaries of the second inequality of IIH are depicted by (o). The 
bright shaded area is obtained from Eq. 11211 . which is the approximation valid in the linear response 
regime. 



the rightward one indexed by ^r, and the leftward one indexed by zl, Eq. © gives 

" ° \{iH.|*iRe(o,to)} {iL\ti-^e{o,to)} J 

where {qi^) = Qr = q'^'^ and {(ji^) = qL = —q°'^ [10]. kr and kl are rate constants for the 
rightward and the leftward stepwise motion respectively and in equilibrium kr = kl = k*"^. 
The expression of Lpq in Eq. (|10() is equivalent to Eq. ®. Furthermore, direct numerical 
calculations of the time integral of Eq. ||SJ) shows a good agreement with q°'^Dp/l (data not 
shown) . 

Qualification of the heat pump. - A system functions as a heat pump if it can absorb 
energy from a cooler heat bath and dissipate it to a hotter one. The present system can 
function as a heat pump when the conditions 

ATJq < 0, I Jql > W/2 > (11) 

are valid. The former inequality means that the heat flow Jq between the two heat baths 
is against temperature gradient AT. The latter means that energy is absorbed from one 
of the heat bath while dissipated to the other, i.e., J1J2 < where Ji = — Jq — W/2 and 
J2 = Jq - w/2. 

For the linear response region with small / and AT, the conditions Hll|) are reduced to 

AT(AT + m'^^f/Xl) < 0, (12) 

which is satisfied in the bright shaded area in f-AT plane of Fig. Here we substituted Jq 
of Eq. lO and Q into the first inequality of l|ll|) . W is vanishing in the linear order of / and 
AT because W = Jpf, which means that the second inequality in ((TT|l always holds. 

For larger values of / and AT, nonlinear effects (e.g. nonvanishing W) are not negligible 
and the region where the system can work as a heat pump is reduced. For instance, for large 
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values of / with AT ~ 0, we obtain Ji < and J2 < as is seen in Fig. |2Ia) although Eq. 
p2|l is satisfied. The region of / and AT satisfying the conditions where the system 
actually functions as a heat pump are numerically identified in the model of Fig. ^a) as the 
dark shaded area in Fig. El Adjusting the strength of the apphed force within this confined 
area of parameters, we can choose the energy absorption rate from a cooler heat bath. The 
rate is maximum at |/| ~ |'z'"'|/^ m the present system (refer Fig. (^). The direction of 
the heat transfer can also be controlled by the direction of the applied force. The rightward 
mechanical force induces the absorption of heat from B2 and the dissipation to Bi (i.e. Ji < 
and J2 > 0), and vice versa for the leftward force (Ji > and T2 ^ 0). 

With this heat pump, we can cool an object to the temperature limited by the conditions 
Hll|l optimizing the applied force. If we need to achieve lower temperature, we should construct 
a proper structure combining multiple pumps. 

Discussion. ~ In this paper, we have shown that heat transfer can be controlled by 
a mechanical force. We clarified the conditions for the mesoscopic system to function as a 
heat pump based on the study of the linear response regime and numerical simulations in fully 
nonequilibrium. These results are consistent with the results obtained in a discretized solvable 
model for Feynman's ratchet [4]. The Langevin dynamics adopted in this study suggests an 
intuitive picture based on the barrier overcoming processes and associated heat transfer q°'^. 
The form of the reciprocity coefficients in Eq. H1U|I generally applies to the class of systems 
with an asymmetric periodic structure and two heat baths. When the direct measurement of 
q^'^ is difficult, we can estimate q'^'^ from the Onsager coefficients as q°'^/l — Lqp/Lpp. This 
estimation can be apphed to other models of heat pump, e.g. [4] and [12]. 

The form of the reciprocity coefficient Lqp = q°'^K°^l implies that the present heat pump 
works effectively in small energy scales where barrier-overcoming processes are governed by 
thermal fluctuations. This is typically realized in molecular scales. The pump is not efficient 
in very low temperatures where the barrier-overcoming event rarely occurs. It would also be 
less effective for very high temperatures. This is because the value of q°'^, which depends on 
the left-right asymmetry of the system, would be smaller in higher temperatures. 
Acknowledgment: We are grateful to H. Tasaki for a critical reading of this manuscript. 
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Note: After the completion of the present work, [13] and [14] were published and [12] 
appeared in the archive, where different models of heat pump are discussed. 
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